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We have developed a method for extracting equilibrium observables from non-equilibrium simulations by rapidly 
changing the temperature and recording the subsequent evolution of the Polyakov loop. Both nucleation and 
spinodal decomposition are observed to occur. In the latter case the Polyakov loop correlation function shows 
exponential growth for wavenumbers less than or equal to the critical wavenumber k c . We have constructed 
the bare as wee as the effective potential for the Polyakov loop, from which k c and mo/kc can be extracted as 
a function of temperature. The shift from spinodal decomposition to nucleation as the dominant equilibration 
mechanism occurs at the spinodal temperature that separates these two regimes. 



1. Introduction 

Spinodal decomposition is characterized by the 
appearance of explosive growth in the low mo- 
mentum modes of the system, while the high mo- 
mentum modes relax to their equilibrium distri- 
butions. The critical wavenumber k c which sep- 
arates low and high momentum behaviors is de- 
termined by the thermodynamic potential which 
governs the long wavelength behavior of the sys- 
tem. We have previously shown that the confined 
phase of pure SU(2) gauge theory decays via spin- 
odal decomposition after a rapid change from be- 
low the deconfinement temperature to tempera- 
tures above [Q. 

The first order character of the deconfinement 
transition in pure SU (3) gauge theory implies 
that the confined phase is metastable for a range 
of temperatures above the deconfinement temper- 
ature Td- At higher temperatures T > T s > Td 
the confined phase is unstable and decays via 
spinodal decomposition. The spinodal temper- 
ature T s separates the metastable and spinodal 
temperature ranges. 

2. Dynamics of the Polyakov Loop 

As a consequence of asymptotic freedom, the 
perturbative form of the effective potential for 
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the Polyakove loop at high temperatures is valid. 
Writing the Polyakov loop in the form Pp = 
1 + 2cos(2ir/3 — ip), the perturbative form of the 
effective action is at 1 loop [|l @] 
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where the constant black body term has been dis- 
carded. 

We assume that the long distance behavior of 
the Polyakov loop 2-point function can be calcu- 
lated using the 1-loop perturbative effective ac- 
tion S e f f and the Langevin equation, 
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where 77, the Gaussian noise term, is such that 
the equipartition theorem holds (l). A linearized 
solution to the Langevin equation gives for the 
Polyakov loop 2-point function (also refered to as 
the structure function) 
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where a = 4TT 2 /g 2 . 

When the equation is linearized about tp = 
the mass squared is negative m 2 = — k 2 = 
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— g 2 T 2 /3 causing exponential growth for values 
of the wavenumber less than the critical value k c . 
At late times, the equation is linearized about the 
non-trivial minimum at tp = 2tt/3. Here the mass 
is the Debye screening mass ran = gT, and the 
system relaxes to the non-trivial minimum. 

3. Method 

3.1. The Bare Potential 

After integrating out all other fields at the scale 
a we are left with an action just in terms of the 
Polyakov loop. The minimal Landau- Ginsberg 
action which is gauge invariant, Z(3) invariant, 
and only involves terms up to the fourth power 
in the Polyakov loop is 

S = J d 3 x{\VP F \ 2 + a 2 (P F P F ) 

+ y(^) 2 + |(P| + (P;) 3 )), (4) 

where P F = tr F (P) is the trace of the Polyakov 
loop in the fundamental representation. The 
Langevin equation without noise is then used as 
an ansatz for the evolution of the low momentum 
modes. In Fourier space this yields an equation 
of the form 

±F[P F ]=c 1 (h?)F[P F ] 

+c 2 (k 2 )T[(P* F P F )P F ] + c 3 (k 2 )F{(P F ) 2 }. (5) 

The behavior of the lattice simulation is fit to 
the above equation for each k mode yielding the 
parameters. A fit of ci(fc 2 ) versus k 2 gives the 
critical wavenumber a 2 = — k 2 . 

3.2. The Effective Potential 

The potential obtained above has a minimum 
which is not the equilibrium value of the vacuum 
expectation value of the Polyakov loop(see figure 
2). This is the bare potential defined at the scale 
a, while the effective potential can be defined as 
the generator at zero momentum of 1PI vertices. 
To find the effective potential, the behavior of the 
1-point function is fit to the ansatz 

V(P cl ) = a 2 \P cl \ 2 + ^-\P cl \ 4 

+ j (P$ + (^) 3 ) (6) 



where P c i = (P). 

Without multiple momentum modes to fit, the 
absolute normalization of the coefficients of the 
effective potential can not be determined. How- 
ever the ratio of the curvature at the nontriv- 
ial minimum to the curvature at zero is equal to 
the ratio of the square of the Debye mass to the 
square of the critical wavenumber 
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which is comparable to the perturbative result for 
this ratio. The minimum of the effective potential 
is located at the final equilibrium value. 

4. Results 

We start from equilibrated lattices of size 64 3 x 
4 at [3 — 5.5 and change the value of (3 to a 
range of values above the deconfinement value. 
We use the Cabbibo-Marinari heat bath algo- 
rithm both for equilibration and for the evolution. 
The large lattice sizes were required to resolve the 
long wavelength modes, and to reduce the effects 
of changing the physical volume. 

Quenches to just above the deconfinement tem- 
perature show evidence of metastability. However 
more extreme changes in the temperature yield 
behavior consistent with spinodal decomposition. 
Typical behavior of the structure function S(k, t) 
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Figure 1. S(k,r) versus Monte Carlo time. 

after a quench to /? = 6.0 for the three lowest 
non-trivial wavenumbers is displayed in figure 1 



for a single run. The initial exponential rise in 
the structure function is evident. At intermediate 
times the growth levels off, with the appearance 
of "shelves" in some of the modes which is pre- 
sumably a consequence of mode-mode coupling. 
At late times the modes relax to their final equi- 
librium values. 




Figure 2. V(P) vs. P 

In figure 2 the bare and the effective potential 
are plotted. As the overall normalization of the 
effective potential cannot be fixed by the present 
methods, the normalization was chosen so that 
both curves could be plotted side by side. These 
potentials were obtained by averaging over sev- 
eral runs and the error bars were estimated using 
a jackknife analysis. Note the significant shift in 
the location of the minimum from the bare to the 
effective potential. This indicates that fluctua- 
tions are large even at temperatures of the order 
of 2T d . 

In figure 3 the value of the critical wavenum- 
ber determined from the bare potential is plotted 
versus temperature for the four temperatures con- 
sidered. The temperature dependence is fit to the 
ansatz k\ = aT 2 + bT d 2 , which is a form suggested 
by a model for the pure gauge thermodynam- 
ics ||. The fit is used to determine the spinodal 
temperature, which is defined by k c (T s ) = 0. The 
spinodal temperature found is T s — 1.29(7)7^ and 
corresponds to (3 S « 5.81. 
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Figure 3. k c (T) vs. T/T d 

5. Conclusion 

Pure SU(3) gauge theory exhibits both 
metastability and spinodal decomposition follow- 
ing a quench to temperatures above the decon- 
fincmcnt temperature. In the spinodal region the 
evolution can be modeled using Langevin dynam- 
ics, and the bare and effective potentials for the 
Polyakov loop field can be constructed. These 
potentials determine k c and mi)/fc c . We find a 
nontrivial temperature dependence for the critical 
wavenumber. As the temperature is lowered the 
critical wave number decreases. The extrapola- 
tion to zero of k c (T) gives a spinodal temperature 
slightly above the deconfinement temperature. 
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